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ABSTRACT 

Globular clusters have linear sizes (tidal radii) which theory tells us are determined by their masses 
and by the gravitational potential of their host galaxy. To explore the relationship between observed 
and expected radii, we utilize the globular cluster population of the Virgo giant M87. Unusually deep, 
high signal-to- noise images of M87 are used to measure the effective and limiting radii of approximately 
2000 globular clusters. To compare with these observations, we simulate a globular cluster population 
that has the same characteristics as the observed M87 cluster population. Placing these simulated 
clusters in the well-studied tidal field of M87, the orbit of each cluster is solved and the theoretical tidal 
radius of each cluster is determined. We compare the predicted relationship between cluster size and 
projected galactocentric distance to observations. We find that for an isotropic distribution of cluster 
velocities, theoretical tidal radii are approximately equal to observed limiting radii for R gc < 10 kpc. 
However, the isotropic simulation predicts a steep increase in cluster size at larger radii, which is not 
observed in large galaxies beyond the Milky Way. To minimize the discrepancy between theory and 
observations, we explore the effects of orbital anisotropy on cluster sizes, and suggest a possible orbital 
anisotropy profile for M87 which yields a better match between theory and observations. Finally, we 
suggest future studies which will establish a stronger link between theoretical tidal radii and observed 
radii. 

Subject headings: galaxies: individual (M87) - galaxies: kinematics and dynamics - globular clusters: 
general 



1. INTRODUCTION 



Historically it has been assumed that the gravitational 
field of the host galaxy re gulates its satellite sizes, includ- 
ing globular clusters (e.g.lvon Hoerner lll957tlKing 1 1 19621 ; 
Innanen. Harris, fc Webbinkl 119831: Uordan et al.l I200E 



Binnev fc Tremamdl2008t iBertin fc Varrill2008D .The size 
of a globular cluster is more commonly referred to as 
the tidal radius, and has both a theoretical and observa- 
tional definition. Theoretically, the tidal radius is better 
known as the Jacobi radius, which marks the distance 
past which, a star will feel a stronger acceleration towards 
the galaxy and escape the globular cluste r. First-order 
tidal theory determines the tidal radius (|von Hoerner I 
1957) via: 



U9c[ 2M > 



(1) 



where R gc is the galactocentric distance of the cluster, 
M is the the cluster's mass, and M g is the mass of the 



For 
to 



galaxy, assumed in early studies to be a point mass 
an isothermal halo, we have M(R gc ) oc R gc and i 

2 

first order we should expect r t cx R gc if the mean cluster 
mass does not vary strongly with galactocentric distance. 
In addition, if the structural properties of the clusters 
such as their central concentrations c do not vary sys- 
tematically with R gc either, then the mean effective (or 

2 

half-mass) radius should also increase as cx R gc . 

A rough illustration can be drawn directly from the 
Milky Way. Taking Milky Way cluster effective radii 
from Harris 1996 (2010 Edition), we plot the mass- 
normalized radius log rh/{M$) versus the log of each 
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cluster's three dimensional galactocentric distance R3D 
(top panel) and projected galactocentric distance R2D 
(bottom panel) in Figure Q] We assume the Y-Z 
plane of the Milky Way i s t he plane of the sky. See 
Ivan den Bergh et al.l (| 1 9 9 If ) or lMackev &: van den Berehl 
( 2005f7 for earlier versions of this plot. The purpose of 
normalizing cluster size by cluster mass is to remove 
any scatter due solely to differences in cluster mass as 
seen in Equation [T] The lines of best fit indicate that 
r h oc ±0 ' 06 and r h oc {R 2D ) OA6±om for a given clus- 
ter mass. Although the projection to 2D does reduce the 
slope of < rh > to something close to i?2, it certainly 
does not flatten the trend entirely. While the relationship 
between cluster and three dimensional galactocentric 
distance is comparable to the prediction from tidal the- 
ory, we do not expect an exact match as the Milky Way 
is not a spherical galaxy, and clusters do not have circu- 
lar orbits as required in the calculation. The remaining 
scatter around the lines of best fit is expected to be due 
at least partly to differences in cluster orbits and central 
concentrations. 

The bottom panel will be what we may expect to see for 
any giant galaxy that is projected onto the plane of the 
sky. In essence, basic tidal theory predicts that clusters 
in the outer halo of a giant galaxy are permitted to have 
much larger linear sizes than those in the inner halo, 
other things being equal. 

However, an intriguing puzzle beginning to emerge 
from recent measurements in a variety of giant E galax- 
ies is that the observed trend is very much shallower 
than expected, near rh oc i?° c 1-0 - 2 , where now we de- 
note R nr as the projected (2D) galactocentri c distance. 
(e.g. iGomez & Woodlevl 120061: lHarris l [2009b). See also 
iSpitler et al.l (|2006f ) and lHarris et aLlpOlOl ) for similar 
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Fig. 1. — Log(r h /M"3 ) versus the log of each cluster's three di- 
mensional galactocentric distance (top panel) and projected galac- 
tocentric distance (bottom panel). The red dashed lines mark the 
linear lines of best fit. 

results from the giant Sa galaxy M104. These results do 
not yet have a full explanation. 

These newly measured trends hint that other factors 
may be in play in addition to the simple tidal theory 
outlined above. In this paper, we explore particularly one 
possible route to explaining the observed shallow trend of 
< rh >, which is to invoke an anisotropy gradient in the 
cluster orbits. // outer halo clusters tend to have more 
strongly radial orbits than the inner-halo ones, they will 
be carried deep into the potential well of their parent 
galaxy and thus trimmed back to smaller-than-expected 
radii. More will be said about this in Section 3 below. 

The observational size of a globular cluster is the limit- 
ing radius, which is "the out er limit of the cluster wher e 
the density drops to zero" (jBinnev fc Tremaind |2008). 
For clusters outside of the Milky Way, the only way of 
defining the exact size of a cluster is to examine its den- 
sity profile, or surface brightness profile. This is done by 
fitting observ ed pro files to models such as the well known 
[Kh7gf (fl96l) rK62 ).lKlngl (fl966h l'K66). IWiison1 ((19751 ) 
rW75). or lSersicI (11968D (S68^ models. 

While it is generally assumed that the Jacobi ra- 
dius and the limiting radius of a cluster are the same, 
recen t studies fe.g iBrosche. Q dcnkirchcn. & G effert I 
19991: iGieles fc BaumgardtJ 120081 : iKupper et all 120101 : 



Baumgardt et al.ll2010f) are finding that this assumption 



may require modification. A way to explore this assump- 
tion is to directly compare the relationship between clus- 
ter size and galactocentric distance for both an observed 
and simulated cluster population. If the correct relation 
can be established, then we can begin to utilize globular 
clusters in new ways to extract further information, such 
as the mass distribution within a galaxy and the orbital 
distribution of the clusters. 

As an initial test case, we use the globular cluster pop- 
ulation around the Virgo giant M87, which provides an 
exceptionally large number of clusters whose sizes are 



well resolved by the Hubble Space Telescope (HST). HST 
Archive images of M87 are available that are unusually 
deep and high signal-to-noise, containing nearly 2000 of 
its globular clusters. From a theoretical point of view 
the gravitational field, which is required for the calcula- 
tion of a cluster's tidal radius, is better k nown for M87 
than for any other giant elliptical galaxy (|McLau ghlin 
1999). In addition, th e foreground redden ing and field 
contamination are low (jTamura et al.ll2006h . Essentially, 
M87 provides the best available testbed for comparing 
theoretical and observational tidal radii. 

In Section 2 we will fit the observed surface brightness 
profiles of M87 globular clusters with different models, 
determine the tidal and effective radii of each cluster, 
and establish the observed trend between cluster size and 
galactocentric distance. In Section 3 we make use of ob- 
servationally determined parameters of M87 to simulate 
a theoretical cluster population. Using the known grav- 
itational field of M87, we can determine each simulated 
cluster's perigalactic distance and calculate theoretical 
tidal radii. Projecting the simulated clusters onto a two 
dimensional plane, we will then finally have a predicted 
relationship between cluster size and projected galacto- 
centric distance. The results of the simulation are then 
compared to the observational results for both isotropic 
and anisotropic distributions. Our conclusions and fu- 
ture work are then discussed in Section 4. 

In what follows we adopt (to - M) = 30.95 for M87, 
to keep consi s tency with the mass distribution model of 
iMcLaughlinl (p99l ). 

2. OBSERVATIONS 

The HST ACS/WFC Archive images we use in this 
study are from program GO-10543 (PI Baltz). The co- 
added composite exposur es in each filter w ere the same 
ones described in detail in lBird et al.l (|2010[). constructed 
throu gh use of the APSIS software (|Blakeslee et al.1 
2003), which performs accurate image registration, 
cosmic-ray rejection, and distortion correction with driz- 
zle. W hile the raw i mages are the same as those 
used bv iMadrid et al.l (|2009f ). iPeng et al.1 (|2009h and 
IWaters et al.1 ( 20091 ) . subpixel resampling was done dur- 
ing the drizzle step to yield final combined science images 
with an improved scale of 0'/025 px _1 (half the native 
pixel size of the camera). As will be seen below, the 
subsampling produced a noticeable improvement in the 
effective spatial reso lution of the data compared with all 
previous studies (see iBird et al.l[2010l for a more detailed 
description) . 

As an initial step, the ELLIPSE and BMODEL func- 
tions within STSDAS were used to fit elliptical isophotes 
to the brightness distribution of M87 in both the V and 
I images. Subtracting the isophotal model from the im- 
ages to remove the brightness of the galaxy allowed for 
the easier detection of objects that were previously hid- 
den (see Ferrarese et al. 2006 or Peng et al. 2009 for 
similar examples). 

The residual images in F814W and F606W were then 
used to identify globular cluster candidates. After de- 
termining the standard deviation of the background sky 
pixel values, we searched for candidates brighter than 
threshold cut-offs chosen to be faint enough to include 
all true clusters but bright enough to exclude almost al l 
individual halo stars in M87 itself. The lBird et~aTl (|2010D 
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Fig. 2. — CMD of the globular cluster candidates in M87. True 
magnitudes were obtained through aperture photometry extrapo- 
lated to large radius and converted to V and I. The familiar blue 
and red sequences are clearly visible. 

analysis discusses the much more challenging measure- 
ment of the faint halo stars. After objects were matched 
based on their position on both images, we manually re- 
moved candidates very near the center of M87 (R gc < 0.2 
kpc or 100 pixels) and near the M87 jet, as the back- 
ground light intensity is much higher in these regions. 
This resulted in 2052 globular cluster candidates. 

To construct a color-magnitude diagram (CMD) of the 
candidates, instrumental magnitudes were converted to 
true magnitudes thr ough aperture photo metry extrapo- 
lated to large radius (jSirianni et al.ll2005l). and converte d 
to (V,I) with the transformations in iSaha et al.l (|201lD . 
The CMD of the candidates is shown in Figure [5J It re- 
veals the familiar blue (metal-poor) and red (metal-rich) 
sequences that are now well established. For a more de- 
tailed ana lysis of th e clust er photometry, we refer the 
reader to lPeng et all (|2009ft . 

With the list of remaining globular cluster candidates, 
the surface brightness distribution of each cluster was 
ht with PSF-convolved K62, K 66, W76, and S68 mod- 
els via the c ode G RIDFIT (e.g. iMcLaughlin et aLl l2008l; 
lHarris et all l2010t iBarmbv et all 120071 ) The code re- 
turned best-fit values of effective radius (r^) and central 
concentration which for the King models is the familiar 
c = log(^r). To remove non-clusters we eliminated ob- 
jects with x 2 values greater than 10 and all candidates 
with c < 0.5 and c > 3.0, which indicates a poor fit to 
the observations. Finally, we removed objects with large 
differences between the effective radius in the V and I 
bands, leaving us with 1290 globular clusters identified 
with high confidence. 

We have carried out size measurements on both the 
original and isophote-subtracted images. Focusing on 
clusters in the outer regions of M87, where the light 
gradient is shallow and the influence of the subtrac- 
tion is minimal, we found the mean offset (r^ (original)- 
rh (subtracted)) to be 0.13 pc. Since the surface bright- 



FlG. 3. — V magnitude distribution of all the globular clus- 
ter candidates (un-shaded) and candidates which are not elimi- 
nated by the x 2 j central concentration and |Ar^(V — I)\ criteria 
(shaded). The usual globular cluster luminosity function with a 
turnover magnitude of M v = — 7.3 and a standard deviation of 1.3 
is shown as a solid line. The luminosity function of the final candi- 
date list is Gaussian-like, symmetric about a visual magnitude of 
M v = —7.6 with a standard deviation of 1.0 (dotted line). 

ness profiles in the subtracted image are internally more 
precise, we take the effective radius as determined by 
model fits to clusters in the subtracted image and add the 
mean offset. This r esults in ou r clust er siz es being com- 
parab le to those of iPeng et al.l ([2009) and lMadrid et all 
(2009), with a mean difference of approximately 0.03 pc 
between both studies. 

To examine which objects were removed from the glob- 
ular cluster candidate list, we consider the V magnitude 
distribution of objects before and after the x 2 , central 
concentration and \Arh(V — I)\ cuts were made. In Fig- 
ure [3] we see that the original dataset (un-shaded) is 
in agreement with the usual globular cluster luminosity 
function, which has a turnover magnit ude of M„ = —7.3 
and a standard deviatio n of 1.3 (e.g. iBrodie fc Straderl 
I2006t IPeng et al.l l2009T i. However, after the candi- 
dates have been removed based on the cuts described 
above, the distribution is reasonably approximated by a 
Gaussian-like function, symmetric about a visual magni- 
tude of M v — —7.6 with a standard deviation of 1.0 (dot- 
ted line) , which we use below to set up our model simula- 
tion. It appears that model fits to the brightness profiles 
of faint clusters provide less accurate central concentra- 
tions or effective radii. We stress that this sub-selected 
luminosity function shape is used only for the purposes 
of setting up our simulated globular cluster distribution. 

Comparing the results of the K62, K66, W75, and S68 
model fitting for the final list of globular cluster candi- 
dates (Figure 0J, we found that for many objects, the 
K66 tidal and effective radii were significantly smaller 
than those predicted by the other models. Inspecting 
the brightness profiles of these specific clusters, we ob- 
served that clusters with brightness profiles that are more 
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Fig. 4. — K62 effective radius vs. K66 effective radius (left panel), 
W75 effective radius (center panel) and S68 effective radius (right 
panel). The dotted line represents equality. 

extended were the source of the discrepancy. This is a 
known fault of K66 models, as they predict a sharp cutoff 
at the tidal radius while W75 models can more success- 
fully fit a more gradual decrease in brightness or density 
near the tidal radius. Since K62 and S68 models are 
not dynamically motivated, they are less affected by how 
cluster brightness behaves near the tidal radius. This 
led us to reject K66 model fits. Next, we inspected the 
mean and root mean square |Ar^(V — 1)\ for each cluster 
as found by all four models. K62 model fits yielded the 
lowest root mean square scatter, suggesting K62 fits have 
the best internal consistency. Therefore we adopted K62 
measurements to compare with our simulations of M87 
cluster sizes. 

The average K62 effective radius of each cluster in the 
V and I bands is plotted against projected galactocentric 
distance in Figure[5J We choose to compare effective radii 
as opposed to tidal radii because effective radii are the 
directly measured quantity where as limiting radii rt are 
only calculated from and c (see lHarris et al.ll2010l for 
further discussion). Looking at the relationship between 
effective radii and projected galactocentric distance, we 
find that the basic trend is much flatter than predicted 
by tidal theory or the one found for Milky Way clusters 
in Figure [1] (bottom panel) . 

3. SIMULATION 

3.1. The Isotropic Case 

We now simulate a globular cluster population orbiting 
within the galactic potential of M87, in order to calcu- 
late theoretical radii as a function of R gc . Each globular 
cluster was given a position in the halo (R, 9, </)), veloc- 
ity (v r , vg, v^), mass (M), and central concentration (c). 
The distribution parameters used in the simulation are 
drawn from Gaussians with parameters summar ized in 
Table [H The spatial distribution was taken from lHarris I 
(2009aj), who found that the projected radial profile of 
the blue and red globular cluster subpopulations could 
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Fig. 5. — K62 effective radius vs. log projected distance for 
observed globular clusters. The solid red line indicates the median 
effective radius calculated using radial bins that are 0.1 kpc in size. 

be fit with a standard Hubble profile relating density 
(<7 C ;) to projected distance (R): 

a cl (R) = <jo/(l + ^r a (2) 

The appropriate values for cto, R , and a are listed in 
Table [T] for the blue and red subsystems. The angular 
distribution was assumed to be spherically symmetric. 

The mass distribution of globular clusters was taken 
from the observed and culled globular cluster luminosity 
function (see Figure [H- Assuming a mass-to- lig ht ra- 
tio of 2 (e.g. iMcLaughlin fc van der Marell (|2005l )). the 
resulting mass distribution is Gaussian about a mean of 
logjjQ = 5.24 and a standard deviation ailogj^) = 0.4. 
In assigning a King-model central concentration param- 
eter (c = log^r) to each globular cluster, the distribution 
of central concentrations in the Milky Way from Harris 
1996 (2010 Edition) was used, a Gaussian distribution 
with a mean of c = 1.5 and standard deviation of 0.4. 

The velocity dispersion (a) in spherical coordinates 
(R, 9, 4>) is taken from the observed line of sight ve- 
locity dispersion of globular clusters in M87 (jCote et al.1 
1200 If ). We have initially assumed the distribution of or- 
bits to be isotropic, such that the anisotropy parameter 

(/3) is zero and an = o-g = a^, where j3 = 1 6 . 2 ^ 

(jBinnev fc Tremaind [2001 ) . For comparis on, the solved 
orbits of Milky Way globular c lusters (jDinescu et al.1 
Il999t ICasetti-Dinescu et al.ll2007t ) do not indicate a pref- 
erence towards circular or radial orbits. In the later dis- 
cussion, however, we relax this assumption and explore 
anisotropic distributions. 

Exactly 6000 globular clusters were simulated. Not 
only does this provide a statistically significant num- 
ber of clusters, but also results in the same number of 
clusters within 10 kpc of M87 as the observed dataset. 
40% were designated as "red" clusters and had positions 
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TABLE 1 

Simulated Globular Cluster Population Input Parameters 



Parameter 



Value 



Radial Distribution 
Blue Population 
co 

Ro 
a 

Red Population 

Ro 
a 

Angular Distribution 
Mass- To-Light Ratio 
Mass Distribution 
(log(M/M )) 

Vlog(M/M ) 
Velocity Dispersion 

(v) 
<j v 
P 

Central Concentration 

(c) 



Hubble Profile 

66 arcmin - 2 
2.0' 
1.8 

150 arcmin -2 
1.2' 
2.1 

Spherically Symmetric 
M/L = 2 
Gaussian 
5.24 
0.40 
Gaussian 
-19 km/s 
401 km/s 


Gaussian 
1.5 
0.4 



drawn from the red Hubble profile from Table [T] The 
remaining 60% were designated as "blue" clusters and 
were drawn from the appropriate observed parameters 
in Table [1] The ratio of total number of blue clusters to 
red clusters is in agreement with the Hubble profiles in 
lHarrisl(l2009l . 

An isotropic distribution of cluster velocities 
has a broad range of orbital eccentricities (e.g. 
Ivan den Bosch et al.1 (fl999)). The normal approach to 
calculating tidal radii is to assume the tidal radius of a 
cluster is imposed at perigalacticon, where the tidal field 
of the host galaxy is th e strongest. This a ssumption 
was in i tially suggested by Ivon Hoerner I ([1957D and later 
iKhTgl ((19621 . and we will use it here. The assumption 
follows from the fact that for almost all real clusters, a 
cluster's relaxation time (t Th ) is greater than its radial 
period, such that a cluster returns to perigalacticon 
before it is abl e to relax. We note, however, that recent 
studies (e.g. iBrosche. Odenkirchen. fc Geffert I 119991 : 
iKupper et all I2010H suggest that instead some sort of 
orbit averaged distance may result in a more accurate 
tidal radius. 

Using the simulated parameters of each cluster, we cal- 
culate t he theoretical tidal rad ius of each cluster as de- 
rived bv lBertin fc Varril (|2008f l: 



n 



Where fL k and v are defined as: 



n 2 = (d<t> G (R)/dR) Rp /R p 



= m 2 + (dH G (R)/dR 2 ) Rp 
v = 4 - k 2 /q 2 



(3) 

(4) 
(5) 
(6) 



orbital frequency of the cluster, k is the epicyclic fre- 
quency of the cluster at R p , and v is a positive dimen- 
sionl ess coefficient . The mass profile of M87 was taken 
from lMcLaughlinl (|1999ft : 



M tota i(r) = M stars (r) + M dark (r) 



M, 



,.(r) = 8.10 x 10 11 M© 



. (r/5.1kpc) 
\l + r/b.lkpcY 



11-67 



(7) 



(8) 



M dark {r) = 7.06xl0 14 M Q x[ln(l+r/560fcpc)- 



(r/560kpc) 



(1 



$G is the galactic potential due to the mass profile of 
M87, R p is the cluster's perigalactic distance, f2 is the 



r/560kpc) 
(9) 

Since clusters are randomly assigned a position, mass, 
central concentration, and velocity, we must also deter- 
mine which clusters are expected to survive to present 
day. More specifically, clusters with orbits that bring 
them to small perigalactic distances may either evaporate 
over a Hubble time due to tidal dissolution or be elimi- 
nated by dynamical friction. Each cluster's tidal dissolu- 
tion time (tdis) was calculated assuming tdi s — 30 x t r h, 
where t r h is a cluster's half- mass relaxation time deter- 
mined with Equation 6 from iMevlan et all ((2001). The 
factor of 30 i s con sistent with NBODY simulations by 
iTrenti et al.1 (|2007h and within the rang e for a , cluster 
in a tidal field from iBinnev fc Tremaind (|2008[) . Clus- 
ters with tidal dissolutions times less than 10 Gyr were 
rejected from the simulation, as they would have fully 
evaporated and not be observable today. Similarly, af- 
ter calc ulating each cluster's infall time due to dynamical 
friction (jBinnev fc Trem ainc 200§), we eliminate all clus- 
ters with infall times less than 10 Gyr. 

To best compare to observations, we first convert all 
tidal radii to effective radii by assuming each simulated 
cluster can be represented by a K62 model. We then 
project the three dimensional position of each cluster 
onto a two dimensional plane, which represents the plane 
of the sky. Each simulated cluster's effective radius is 
plotted against its projected galactocentric distance in 
Figure [S] For comparison purposes, we have also plotted 
each observed cluster in red and the observed median as 
a solid black line. 

3.2. Anisotropic Cases 

Up to this point, we have been operating under the 
assumption that the velocity distribution of M87 is 
isotropic (0 = 0). However, the value of may differ 
in diff erent regions of a galaxy. For M87, ICote et aLl 
(|2001l) found that while the cluster population of M87 
appears to be isotropic as a whole, the inner regions of 
M87 could possibly have a negative value for such that 
orbits are preferentially tangential (0 < 0), while the 
outer regions of M87 m ay contain cluste r s wit h more 
radial orbits (0 > 0). (Weijmans et all (|2009l ) made 
similar conclusions from observations of NGC 3379 and 
N GC 821. Theoreti c al wo rk regarding the Milky Way 
by [Pri cto & Gncdin (2008) an d on d ark matter haloes 
bvlZait. Hoffman, fc Shlosmanl (|2008l ) and iLudlow et al.l 
(2010) among others all find that increases from the 
isotropic case in the inner regions of a galaxy to more 
radial orbits in the outer regions. Additionally, many of 
these same studies all find evidence for to be less than 
zero at small galactocentric distances. 
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Fig. 6. — Effective radius of each simulated globular cluster 
(blue) compared with its projected distance. The dashed black 
line indicates the median effective radius calculated using radial 
bins that increase in size at a rate of 10 0,4 kpc. The observed 
clusters (red) and median (solid black line) from Figure \5\ are also 
plotted. 

The simulation outlined above was repeated for dif- 
ferent values of /?. For (3 greater than zero, the radial 
velocity dispersion was ass umed to be equa l to the ob- 
served velocity dispersion (|C6te et al.ll200l[) . and the 9 
and <f> distributions were assumed to be equal (ag = o^) 
and calculated via the /3 equation. Since o~r > ag, clus- 
ters are brought deeper into the tidal field of the galaxy 
to smaller perigalactic distances. This in turn reduces 
their tidal and effective radii. For (3 less than zero, the 
8 and <j) velocity distributions are assumed to be equal 
to the observed velocity distribution, with the radial ve- 
locity distribution determined via the (3 equation. Since 
o~e > o~r, clusters will have large tangential velocities and 
small radial velocities, keeping them far from the galactic 
center. 

4. COMPARING THEORY AND OBSERVATIONS 

4.1. The Isotropic Case 

The first comparison that was made between observa- 
tional and theoretical cluster radii was for the simulated 
isotropic case, with (3 = 0. The median effective radius 
versus projected galactocentric distance is shown as the 
dashed black line in Figure [5] Simulated clusters (blue) , 
as well as the observed clusters (red) and observed me- 
dian line (solid black line) are also plotted. The first ob- 
servation that can be made is that the simulation and ob- 
servations agree in the inner regions of M87. The second 
observation which can be made is that in the mid-halo to 
outer regions of M87, the simulation appears to overes- 
timate observed cluster sizes. However, at larger galac- 
tocentric distances theoretical tidal radii may be larger 
than observed tidal radii, as it is possible that some clus- 
ters may be tidally under-filling. If they formed at a large 
galactocentric distances where the tidal field is weak, and 
it is likely that they will stay tidally under-filling for the 
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Fig. 7. — Effective radius of each globular cluster compared with 
its projected distance for different values of f}. The dotted black 
line is the median effective radius taken from K62 model fits to the 
observations (Figure l5t . 

duration of their lifetimes (jGieles et al.ll2010l) . Due to 
the limited field of view of the observations, we cannot 
yet make comparisons past R gc ~ 9 kpc. 

4.2. Anisotropic Cases 

We next compare the observational results to simula- 
tions with — 1 < (3 < 1 as shown in Figure [7] with the 
median observational K62 effective radii shown as a dot- 
ted black line. Along each of the model curves, (3 is held 
constant with R gc at the value shown (that is, there is 
no /3-gradient). Simulations with (3 = —0.5 and f3 = 0.2 
were also performed. As expected, the /? = —0.5 results 
are between the j3 = — 1 and j3 = median lines in Fig- 
ure El while the j3 = 0.2 results are between (3 = and 
(3 = 0.5. 

Since none of the single (3 simulations result in a per- 
fect match between theory and observations, Figure [7] 
suggests that an anisotropy profile j3(R gc ) must be in- 
corporated into the simulation. Matching simulated and 
observed effective radii suggests that the globular cluster 
population is approximately isotropic ((3 = 0) for R gc < 2 
kpc, and radially anisotropic ((3 > 0) for R gc > 2kpc. To 
determine a possible anisotropy profile for M87, we per- 
form a x 2 test between the observations and simulation 
in different radial bins. 

We combine our findings in Figure[7]with the % 2 testing 
to create a possible anisotropy profile of M87 such that 
13 = for R gc < 2 kpc and /3 = 0.8 for R gc ~ 10 kpc. We 

chose to simulate cluster populations with (3 increasing 

i 

proportional to R gc beyond 2 kpc. The effective radius 
distribution produced by this profile is plotted in Figure 
[H For convenience, we have also plotted the results of 
the (3 = simulation from Figure [5] in cyan. While this 
possible anisotropy profile yields a stronger agreement 
between theoretical and observational tidal radii, we view 
this profile as only a preliminary one and a broader range 
of anisotropy profiles should be explored. Nevertheless, 
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Model: (3 = 
Model: (i a R 
Observations 




Log (R go ) (kpc) 

Fig. 8. — Effective radius of each globular cluster compared with 
its projected distance for a radially dependent /3(r) (red) and for 
/3 = (cyan). The black line is the median effective radius taken 
from K62 model fits to the observations 

our exploration of the range of /3-models indicates that 
the imposition of a /3-gradient is capable of matching the 
observed r^(i? gc ) trend seen so far. 

Our findings that j3 may increase with galactocentric 
distance are generally consistent with previous observa- 
tional and theoretical results. However, having (3 in- 
crease very quickly to a value of 0.8 at 10 kpc has not 
been found by previous studies. With a more gradual 
change in /?, we would see that theoretical Jacobi radii 
are greater than observed K62 tidal radii in the mid- 
halo region. This can partially be attributed to clusters 
under-filling their tidal radius, as discussed in Section 
4.1. In the inner regions of M87, where the tidal field is 
much stronger and globular clusters are expected to be 
tidally filled, we see a strong agreement between obser- 
vationally determined and theoretical calculated cluster 



5. CONCLUSIONS AND FUTURE WORK 

We have measured accurate effective radii of approxi- 
mately 2000 globular clusters within 10 kpc of the center 



of the giant elliptical galaxy M87. A theoretical clus- 
ter population was then simulated with the same ob- 
servational characteristics of M87 and the tidal radius 
of each cluster was d etermined using the formalism of 
iBertin fc Varrl (|2008l) . 

The relationship between median cluster size and pro- 
jected galactocentric distance was used to compare theo- 
retical and observational tidal radii. To first order, it 
appears that the observational and theoretical cluster 
sizes are the same. Upon closer inspection, the theo- 
retical and observational distributions are not in com- 
plete agreement, as tidal theory tends to overestimate 
cluster sizes in the outer regions of M87. Unfortunately, 
comparisons in the outer regions of M87 are constrained 
by the radial limit of our observations. In an upcoming 
HST Cycle 19 program, we will be able to add cluster 
size measurements extending to R gc ~ 80 kpc, where 
stronger tests of the theory can be made. 

One possible explanation for the discrepancy between 
theory and observations that we explore here is the ef- 
fect of orbital anisotropy on the simulated distribution 
of cluster sizes. We compare simulations with different 
values of (3 to observations, and conclude that M87 may 
have an anisotropy profile with /3 ~ in the inner regions 
of the galaxy and j3 > in the outer regions. This sort 
of profile is consistent with recent observational and the- 
oretical results. We include a possible anisotropy profile 
which yields a stronger agreement between theoretical 
and observational tidal radii. 

Some issues about the cluster orbits and their inter- 
nal dynamic evolution clearly remain to be investigated. 
Future work will include the use of Monte Carlo Markov 
Chain formalism to explore a broad range of anisotropy 
profiles to find the best possible match. Future N-body 
work will also explore the assumption that a cluster's 
tidal radius is imposed at perigalacticon. 
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